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ABSTRACT 

Glitches are sudden changes in rotation frequency and spin-down rate, observed from 
pulsars of all ages. Standard glitches are characterized by a positive step in angular 
velocity (Afl > 0) and a negative step in the spin-down rate (All < 0) of the pulsar. 
There are no glitch-associated changes in the electromagnetic signature of rotation- 
powered pulsars in all cases so far. For the first time, in the last glitch of PSR J1119- 
6127, there is clear evidence for changing emission properties coincident with the glitch. 
This glitch is also unusual in its signature. Further, the absolute value of the spin-down 
rate actually decreases in the long term. This is in contrast to usual glitch behaviour. 
In this paper we extend the vortex creep model in order to take into account these 
peculiarities. We propose that a starquake with crustal plate movement towards the 
rotational poles of the star induces inward vortex motion which causes the unusual 
glitch signature. The component of the magnetic field perpendicular to the rotation 
axis will decrease, giving rise to a permanent change in the pulsar external torque. 


1 INTRODUCTION 

The 2007 glitch of PSR J1119-6127 is unusual and inter¬ 
esting as the first case with clear indications of changing 
pulsar emission properties coincident with the glitch. The 
event is also unusual in its long term signature of decreased 
spin-down rate. These signatures require an extension of the 
vortex creep model which has become the standard model 
for evaluating glitches and post-glitch response. The exten¬ 
sion of the model must also make allowance for changes in 
the pulsar torque suggested by the glitch related changes in 
emission properties. 

Glitches are sudden increases in the rotation rate of 
pulsars followed by relaxation towards the pre-glitch state. 
The fractional change of the angular velocity, Afl/fl, in a 
glitch is in the range ~ 10“^'^ — 10~®. Glitches are usually 
accompanied by jumps in the spin-down rate, Afl/fl, in the 
range ~ 10“^ — 10~^. To date, about 400 gli tches have been 
observed in more th an a hundred pulsars llEspinoza et al.l 
|2011 | : IYu et al.ir2013l ). Since the earliest glitch observations, 
the sudden changes in rotation frequency and spin-down 
rate were seen to relax back towards the pre-glitch values 
on timescales of days to years. This is interpreted as a sig- 
nature of superfluid interior co mponents of the neutron star 
llBavm. Pethick fc Pineal 19691 ). as a star composed of nor¬ 
mal matter would relax much faster. 

Several models have been proposed to explain glitches 
and post-glitch r elaxation. In the early starquake model 
llRudermanI 1 1969f ). the solid crust of the neutron star oc¬ 
casionally cracks under stresses induced by the ongoing 
spin-down of the star, thereby readjusting to a less oblate 
shape closer to the equilibrium shape that a fluid star 


would follow while spinning down. By conservation of an¬ 
gular momentum, the reduction in moment of inertia of the 
crust is accompanied by an in crease in its angular velocity . 
Glitches in the Grab pulsar ( Wong, Backer fc Lvnd 1200 ll ) 
and PSR J0537-6910 (iMiddleditch et al. 20061 ) can be ex¬ 
plained by this model. However, starquakes cannot explain 
large glitches th at repeat every few years, as exhibited by 
the Vela pulsar dBavm fc Pinesl[l97ll ). The required rate of 
dissipation of elastic energy stored in the solid crust would 
produc e an X-ray luminosi ty enhancement which is not ob¬ 
served (iGiirkan et al.|[2000l ). 

The standard model for the pulsar glitches is the 
vortex pinning—unpinning (vortex creep) model based on 
the dynamics of the neutron star’s su perfluid interior 
dAnderson fc Itohlll975l : lAlnar et al.lll984l ). This model in¬ 
vokes the minimal storage and dissipation of energy for a 
star with angular momentum. The expected energy dissipa¬ 
tion in a large glitch, at the expense of the rotational kinetic 
energies of the two components, does not violate any obser¬ 
vational upper bounds. Models based on pinned superfluid 
components can explai n the various modes of glitch and 
post-glitch behaviour (Haskell. Pizzochero fc Sidervi I2OI2I : 


iHaskell fc Antonopoulou 12014). 


Radio pulsar glitches observed up to the 2007 gl i tch o f 
PSR J1119-6127 IWeltevrede, Johnston fc Espinod I2OIII ) 
showed no glitch correlated changes in the electromagnetic 
signatures, like pulse shape, emission pattern, spectrum and 
polarization. Previous applications of the vortex creep and 
starquake models assumed that there were no changes in 
the pulsar torque at the time of the glitch. Glitches and 
post-glitch response were explained entirely in terms of the 
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internal structure and dynamics of the neutron star. The 
2007 glitch of PSR J1119-6127 shows clear evidence for 
changing emission properties induc ed by the glitch, sw itch- 
ing on intermittent pulses (see, e.g.. lKramer et al]l2006ll and 
also showi ng rotating radio transient (RRAT) behaviour 
(see, e.g., iKeane fc McLauehlinl 1201 il l. Interestingly, this 
glitch also displayed AQ > 0 after transients have de¬ 
cayed, in contrast to the signatures of “standard” glitches 
which are characterized by a negative step in spin-down 
rate (AQ < 0). The high magnetic field radio pulsar PSR 
J1846-0258 had comparable glitch-in duced emission changes 
llLivingstone. Kaspi fc Gavriilll2oToll . The radio pulsar PSR 
J0742-2822 showed a suggestive connection between chang¬ 
ing radio emission and pulse shape features and glitch activ¬ 
ity, however there is currently little direct evidence to estab¬ 
lish a robust link between the m due to absence of enough 
data following the glitch date JKeith, Shannon fc JohnstonI 
l2013h . The RRAT J1819-1458 was also reported to have an 
incre ase in its activity associated with a glitch (iLvne et al.l 

In this paper we analyze the 2007 glitch of PSR J1119- 
6127 and extend the vortex creep model to include the pos¬ 
sibility of a sudden change in the pulsar torque associated 
with the glitch, as suggested by the changing emission prop¬ 
erties, and to address the atypical glitch signature. In (J2]we 
summarize the unique properties of the 2007 glitch of PSR 
J1119-6127. In !j3]we review the vortex creep model, while in 
^we develop the model to include the unusual signatures in 
the spin frequency and spin-down rate, and allow for glitch 
associated changes in the pulsar torque. We apply our ex¬ 
tended model to the peculiar glitch of PSR J1119-6127 in 
j}51 We discuss our results in 


2 THE PECULIAR GLITCH OF PSR 
J1119-6127 

PSR J1119-6127 is a young pulsar with a period P = 
0.41 s and a period deriva t ive P = 4 x 10“^^ Hz s“^ 
discovered by ICamilo et al.l 1 I 2 OOOI I. It has a characteris¬ 
tic age Tc = P/{2P) = 1625 years, and a high surface 
dipole magnetic field B 8.2 x 10^^ G (at the poles) . 
This pulsar has exhibited three glitches llCamilo et al.l[2000l : 

IWeltevrede. Johnston fc Espino^ |2011^ . The third glitch, 

which occurred in 2007, was quite unusual in a number of 
ways dWeltevrede. Johnston fc Espinozall201ll l: 

(i) For a while after the initial exponential relaxation 
is completed, the pulsar is found to be rotating with a 
smaller angular velocity as compared to the pre-glitch value, 
Afl(t) < 0. In the latest data A^l(t) > 0, and may be set¬ 
tling at a positive value IIAntonoDouIou et al.ll2014l l. 

(ii) In the long term, the pulsar slows down at a lower 
rate; the absolute value of the spin-down rate is less (the 
frequency derivative is greater) than its pre-glitch value. All 
> 0 . 

(iii) While the fractional changes in angular velocity are 

small, of the order of 10~^, for the glitches of the Crab pul¬ 
sar, Vela and older pulsars undergo large glitches of size 
Afl/fl as well as smaller “Crab-like” events. The 

2007 glitch of PSR J1119-6127 is a “Vela-like” giant glitch 
from a young pulsar comparable to the Crab pulsar in char¬ 
acteristic age. 


(iv) The radio emission properties of PSR J1119-6127 
displayed changes associated with the glitch. The pulsar 
switched on intermittent pulses and also showed RRAT be¬ 
haviour which seems to have emerged with the glitch. This 
anomalous emission behaviour of PSR J1119-6127 was ob¬ 
served for about three months following the 2007 glitch. 

The mnch smaller second glitch which occnrred in 
2004 may have had similar signatures in the long term 

J iost-glitch frequency and freqnency derivative remnants 
Antonopoulou et al.l 12014 ). The data is sparse, and post- 
glitch evolntion may have been interrupted by the arrival of 
the 2007 glitch. Furthermore, no glitch associated changes in 
emission properties were observed for the 2004 glitch. Here 
we address only the 2007 glitch. 


OVERVIEW OF THE VORTEX CREEP 
MODEL 


The vortex cree p _ model dAlpar et al.l 1 19841 : 

lAlpar. Cheng fc Pined Il989l l attempts to explain the 
processes which cause both the glitches and the post-glitch 
relaxation in terms of a number of distinct superfiuid 
regions in the inner crnst. The superfinid core of the 
neutron star is conpled to the external torque on very short 
timescales, via electron scat tering off magnetized vortices 
dAlpar. Langer fc Saulsiri984ll . The core superfinid therefore 
behaves as part of the effective normal matter crnst. Hence 
the snperfluid component relevant for glitch and postglitch 
dynamics is the crust superhuid. A description of the core 
superfluid blue as well as the crustal superfiuid in terms of 
mutual friction forces acting upon vortex lines is given by 
[Andersson, Siderv fc Comer I (I 2 OO 6 I '). 

The dynamics of the crnst superfiuid is constrained 
by the pinning of the quantized vortex lines to nn- 
clei, interstitial positi ons and po s sibly other strnctures 
in the crust lattice (lAl2ml_[l97^jLink&Epstein|J[l99l|; 
Mochizuki et al. 19991: Avogadro et al. 20081 : Pixzocherc 


20 ^; iHaskell, Pizzochero fc SidervI l2012l : ISeveso et al 


2 OI 4 II . When vortices pin to nuclei, they move with the 


crust’s velocity. A lag a; = fls — flc builds up between the 
superfiuid and the crustal angular velocities fls and flc as 
the crust spins down under the external pulsar torque. This 
lag is sustained by the pinning forces acting upon the vortex 
line. In the case of rotational (cylindrical) symmetry, the 
magnitude of the required pinning force (per unit length) is 
/ = rpsKLO = rpsK{^ls — He), where r is the distance from 
rotation axis, ps is superfiuid density, k is the quantum 
of vorticity. The critical (maximum) lag, uicr, determined 
by the maximum available pinning force, is given by 
oJer = Ep/rbpsK^. Here Ep is pinning energy, ^ is the vortex 
core radius and b is the distance between successive pinning 
sites along the vortex line. If local fluctuations in vortex 
density and superfiuid velocity raise uj above ujcr, there 
will be sudden unpinning and outward motion which can 
lead to an avalanche of vortex discharge (lAnderson fc ItolTI 
I 1975 II . By conservation of angular momentum this leads to 
speeding up of the crust, AHc > 0, observed as a glitch. 
The possibility of such vortex unpinning avalanches taking 
place spontaneously was conh r med by computer simulations 
jMelatos_&_^^szaws^ 20091: Warszawski fc Melatosll201ll : 


IWarszawski. Melatos fc Berlofil 20121 ). 
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Apart from the discontinuous angular momentum im¬ 
parted to the crust by sudden vortex unpinning at glitches, 
the superfluid also spins down continuously between glitches 
by outward flow of vortices. The crustal neutron superfluid 
follows the spin-down of the crust by means of thermally 
activated outwar d creep of vortex l i nes against the pinning 
energ y barriers dAlpar et al.l 1 19841 : lAlpar. Cheng fc Pines! 
1198911 . 

In terms of a simple two component model, involving 
the crust and the superfluid component, the observed spin- 
down of a neutron star’s crust satisfies the equation. 


Ictlc — A^ext + Nint — A^ea;i “ Is^s, (1) 

where Next = /floo is the external torque on the neutron 
star which tries to slow down the crust, and Nint is the 
internal torque arising from the coupling of the superfluid 
to the crust by vortex creep and tends to speed up the crust. 
7c is the moment of inertia of the effective crust (including 
the superfluid core of the neutron star). Is is the moment 
of inertia of the pinned superfluid, while their spin-down 
rates are Cle and Cls, respectively. The spin-d own rate Cls 
of th e superfluid is determined by vortex creep dAlpar et al.l 
[1983). The system reaches a steady state when both the 
superfluid and the crust spin-down at the same rate floo = 
Next/{Is + 7c), sustained at the steady state lag ujao- 

Glitches set the system off from steady state. Post-glitch 
relaxation is due to the recovery of vortex creep, as the lag ui 
inevitably builds back towards steady state due to the ongo¬ 
ing spin-down of the crust under the external pulsar torque. 
The internal torque is so sensitively dependent on the pin¬ 
ning energy Ep and the crustal temperature T that we ex¬ 
pect vortex lines in the different regions of the superfluid to 
respond differently. Depending on the temperature and the 
local pinning parameters in relation to the external torque, 
vortex creep can have a linear or no nlinear dependence on 
the lag dAlpar. Cheng &: Pined [i989l 'n]. In the linear regime, 
the response is linear in the glitch-induced perturbation to 
the lag uj and gives simple exponential relaxation. The relax¬ 
ation time Ti is very sensitively dependent on Ep/kT, with 
Ti oc exp(Ep/kT). The steady state lag cdao = IDIooT; is al¬ 
ways much less than ujcr in this regime. From glitch observa¬ 
tions, up to four ex ponential relaxation terms are seen from 
a particular pulsar dPodson. Lewis fc McCullochll2007ll . 

In the opposite regime we have a very nonlinear re¬ 
sponse to perturbations. The response of a non linear creep 
region k to the glitch will be dAlpar et al.lll984 1. 






1 

1 -f {e^O,k/'r^l 


l)e-‘/Tni 


( 2 ) 


At the time of glitch, creep in those regions which show 
nonlinear response can stop temporarily. These regions de¬ 
couple from rest of the star, so that external torque acts on 
less moment of inertia. Creep restarts after a waiting time of 
to = (5cu/|f2|oo, since the external torque restores the glitch 


^ The claim that the linear regime of vortex creep is never re¬ 
alized for realistic pinning parameters dLin 31201 3) depends on 
the velocity of unpinned vortices, relying on the assumption that 
they move with the global averaged superfluid velocity with drag 
forces, and are not affected by the contributions of interactions 
with the adjacent pinning sites to the local superfluid velocity. 
This issue will be addressed in a separate work. 
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induced decrease in angular velocity lag. The relaxation time 
is 

kT UJer 

In those superfluid regions through which the avalanche 
of vortices unpinned at the glitch pass, moving rapidly 
in the radially outward direction, the ensuing reduction 
5Qs in the superfluid rotation rate determines the offset 
Sljj = dfls -I- AQ,e in the lag, as SQs Aflc This results 
in the response given in Eq. (2), characterized by the wait¬ 
ing time to = 5S7s/|fI|oo > Tni. There can also be nonlinear 
creep regions through which no unpinned vortices pass at 
the glitch, so that Sut = AIIc. In this case to — AS7c/|D|oo 
can be much shorter than r„;, and the contribution of such 
a nonline ar creep region reduces to s imple exponential re¬ 
laxation JCugercinoglu fc Alpaill2M4l 'l. 

= (3) 

like in the case of linear creep regions, but with the nonlinear 
creep relaxation time r„i. 

If we integrate Eq. (2) with the assumption that the 
post-glitch superfluid angular velocity decreases linearly in 
r over the region, correspond ing to uniform de nsity of un¬ 
pinning vortices, one obtains (lAloar et al.ll 1983) 


Ailejt) 7^ r 1 In [l + - l)e ] ] 

^1^ i_e“^ J 

(4) 

In the limit to ^ Tni this reduces to recovery with a constant 
lie 


as observed in the Vela pulsar dAlpar et al.l Il993ll and in 
most Vela-like giant glitches in older pulsars (IYu et al.l 
l2013ll . In the above equations to is the maximum waiting 
time, IA is the moment of inertia of the vortex creep region 
A where unpinning of the vortices has taken place during 
the glitch. Vortices unpinned in regions A pass through re¬ 
gions B with moment of inertia Ib before repinning in an¬ 
other creep region A. Regions B do not participate in spin- 
down by creep, as they do not sustain pinned vortices. Re¬ 
gions B contribute to the angular momentum transfer only 
at glitches, when an avalanche of unpinned vortices moves 
through them. These regions A and B determin e the glitch, 
interglitch and long term behaviour of pulsars dAlnar et al.l 
1199.1 1 199(dl . 

After the exponential transients are removed, observ¬ 
able variables associated with glitches are related to the 
model parameters by the following simple three equations 
dAlpar fc Bavkalll2006l h 


IcAQc = {Ia/^ + Ib)SQs- 

(6) 

ACIc I a 

(7) 

lie I ■ 

Qe = 

I SQs 

(8) 
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Eq. (6) simply states angular momentum conservation and 
gives the glitch magnitude. This is proportional to the num¬ 
ber of vortices which participated in the glitch event. For a 
uniform array of vortices the number of unpinned vortices 
moving outward through radius r is related to the change in 
angular velocity of the superfiuid at r, 

SN = 2TTr^5fls/n — 2nR^5fls/K, (9) 

since r = R, the radius of the star, in the crust superfluid. 
The angular momentum transfer depends on SQa and the 
moment of inertia of the regions that vortices pass through, 
IA and Ib- Eq. (7) is about the torques acting on the pulsar. 
Before the glitch, in steady state, the crust superfluid and 
the rest of the star spin down at the same rate. When a 
glitch occurs, some part of the crustal superfiuid decouples 
from the external torque leading to a jump in spin-down 
rate. Solving these equations for the three unknowns, 1a, Ib, 
and SHa, one can obtain model parameters uniquely without 
making any further assumptions. 


4 EXTENSION OF THE VORTEX CREEP 
MODEL 


In the standard vortex unpinning-creep model only the out¬ 
ward motion of vortices is considered. This gives a nega¬ 
tive post-glitch offset (an increase in the absolute value) 
of the spin-down rate from its pre-glitch value, ACI < 0. 
The spin-down rate relaxes back to the pre-glitch value 
(Afl —>■ 0) for all modes of vortex creep which supply the 
internal torques from the superfiuid acting on the normal 
matter crust. Thus, (i) glitches with the “wrong” sign in 
frequency and spin-down rate require inward vortex motion 
at the glitch; and (ii) long term (persistent) shifts in the 
spin-down rate require either a structural chang e in the neu¬ 
tron star crust, as proposed for the Crab pulsar dAlpar et al.l 
1 19961). or a glitch associated shift in the external torque 
I Link. Epstein fc Bavm|ll992l ). 

Occasional inward fluctuations of vortices, facing an ex¬ 
tra potential barrier, is a low probability component of the 
creep process. Therefore, bulk spontaneous inward motion of 
an avalanche of unpinned vortices is thermodynamically im¬ 
possible in an isolated superfiuid. Large numbers of vortices 
could be transported inward only if the glitch were induced 
by an agent external to the superfiuid, like a starquake. 

Inward vortex motion will increase the superfiuid ve¬ 
locity by some SU'a in regions of superfiuid through which 
vortices have moved inward. Its effect can be investigated 
by changing to with —t'o, where t'o = 5ns/|fI|oo. With this 
we obtain: 


AQa = -^inic 


1 - 


1 -b - l)e-‘/<i 


( 10 ) 


where the primes indicate parameters associated with in¬ 
ward vortex motion. This equation describes the response 
to inward motion of unpinned vortices. When vortices travel 
inward, superfiuid rotates faster. The lag lj thereby increases 
from its steady state value, and creep will be more efficient 
than in steady state, with an enhanced vortex current in the 
radially outward direction. If we integrate Eq. (10) over a 
nonlinear creep region throughout which a uniform average 


density of vortex lines unpinned, or repinned, we obtain: 


Atlcit) 

nc 



1 + (^ni/io) In 


1 -b 


t 



t 



( 11 ) 

The internal torque contribution given in Eqs. (10) and 
(11) leads to an initial positive contribution to Ailc, which 
asymptotically decays to zero. Unlike the nonlinear creep re¬ 
sponse to glitch associated outward vortex motion, as given 
in Eq. (2), the nonlinear creep response to inward vortex 
motion, does not have a waiting time. Instead Eqs. (10) and 
(11) display quasi-exponential relaxation. A constant second 
derivative Clc is not obtained from Eq. (11) when t'o ^ r^i or 
in any other limit. As the integrated response in Eq.(ll) is 
very similar, Eq. (10) is adequate to describe the spindown 
rate when vortices have moved inward. 

Allowing for the starquake induced inward vortex mo¬ 
tion at the glitch, in addition to the natural outward motion 
of many unpinned vortices, we get the following equation in¬ 
stead of Eq. (6), 


IcAnAO) = ( 7 a / + lB)Sna - {IaG + iB'Wa- ( 12 ) 


where / = 1/2 is for the integrated response, Eqs. (4), (5) 
and (11), and / = 1 for the simpler response, Eqs.(2) and 
(10). The hrst term on the right hand side is the angular mo¬ 
mentum transfer due to outward moving vortices, while the 
second term is the contribution of inward moving vortices. 
The physical meanings of Ia' and Ib' are similar to their 
non-primed counterparts. A plate of the crustal solid that 
moves inward in a quake could carry vortices with it, in the 
inward, —r, direction. Nonlinear creep regions with moment 
of inertia I a', and vortex free regions with moment of iner¬ 
tia Ib' are at radial positions between the original and the 
new positions of the plate, and therefore experience a sud¬ 
den increase Sit'a > 0. As creep relaxes back to steady state, 
the net angular momentum transfer from the regions A and 
A! is zero, while the regions B and B' transport angular 
momentum only at glitches and will contribute a remnant 
frequency offset Aitp-. 

IcAflp = 7s5f2s — iB'Sfl'a. (13) 


Extending Eq. (7) to describe the net glitch in the spin- 
down rate with the terms of opposite signs describing the 
response of creep to outward and inward vortex motion, we 
obtain 

Aflc _ Ia_ _ 7^ ,.. .. 

Cla ^ I I ' 

For PSR J1119-6127 the post-glitch Af2 > 0 persists for 
~ 2500 days, as far as the pulsar h as been observed since the 
glitch (lAntonopoulou et al.l 120141) . Here we pursue the as¬ 
sumption that ACIc > 0 is permanent; that it will not decay 
on long timescales in the future. This is a viable assumption 
with the present data, as discussed in the next section. With 
this assumption the permanent shift AClp could be due to 
a structural change in the star, as postulated for the per¬ 
sistent shifts in spin-d own rate observed to accompany the 
Crab pulsar glitches dAlpar et al.lll996l ). or, alternatively, 
due to a glitch associated permanent change in the exter¬ 
nal torque. Unlike the Crab pulsar, the 2007 glitch of PSR 
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J1119-6127 has strong indications that actually the external 
torque has changed, since the pulsar has switched to inter¬ 
mittent and RRAT behaviour with the glitch. It is likely that 
structural changes experienced by PSR J1119-6127 lead to 
a permanent change in the external torque. 


5 MODEL FITS 


We apply a model which is an extension of e arlier appli¬ 
catio ns of the v ortex creep mode l to the Vela jAlpar et al.l 
I1993I I and Crab (lAlpar et al.lll99^ 1 pulsars’ glitches. We take 
one nonlinear creep region with relaxation time T 2 corre¬ 
sponding to the outward motion of the glitches. The new 
component in the extended model is the inclusion of inward 
moving vortices in the glitch, which move through a nonlin¬ 
ear creep region with relaxation time ri (cf. Eq. (10)). We 
also employ a region in which relaxation occurs exponen¬ 
tially with a timescale rs, discussed below. Finally, we in¬ 
clude a possible external torque change as a constant offset 
to the spin-down. We tried model fits with the integrated 
response, Eqs.(4) and (11) and with the simple response 
Eq.(2) and (10). As the residuals are comparable, we choose 
to employ the simple model. 

The expression used for the fit including our extended 
formula is: 


AClc{t) = — ai 


- 02 


1 - 

1 - 


1 -I- aie-0+'^)M 
1 


1 -I- a2e“0+^)/T2 


- + b. 


(15) 


The various parameters are defined by: ai = -^ISIIck), 02 = 
m = (e-*o/"^ - 1), 02 = (6^“/"= - 1), 03 = 
and b = {ANext/N)Cloo, and t is the time since the first 
post-glitch observation, with the time lag A between the 
actual glitch date and the first post-glitch observation. We 
have 9 free parameters. Parameters with the subscript “1” 
denote the contribution from the response of vortex creep 
to glitch associated inward vortex motion, while those with 
subscript “2” and “3” are associated with creep response to 
glitch associated outward vortex motion. 

The exponential relaxation term with amplitude 03 
might describe the response of either an intrinsically lin¬ 
ear creep region, or a nonlinear creep region where there 
was no vortex motion at the glitch, so that the angular ve¬ 
locity of the superfiuid remains unchanged and the glitch 
induced pe rturbation to the angular ve locity lag is simply 
Su = Aflc dGiigercinoglu fc Alparj[2014l ). We adopt the lat¬ 
ter interpretation. This assumption is consistent with the 
results obtained from the fits. 

The moments of inertia of nonlinear creep regions con¬ 
tributing to the long term response are obtained from the 
fit parameters ai and a 2 . The terms a\ and 02 yield the 
numbers of vortices moving inwards and outwards, respec¬ 
tively, during the glitch, b is the long term offset of ACIc 
after all the contributions from creep regions relax back to 
zero. We interpret this as the contribution of the change 
in the external torque. The terms with subscripts “2” and 
“3” contribute AClc{t) < 0 while the parameter b (exter¬ 
nal torque change) and the term with subscript “1” (inward 


motion of vortices) contribute AQc{t) > 0 (see Figure 1). 
This term has the longest time constant, n ^ T 2 > t^. The 
data could also be fitted by assuming no change in the ex¬ 
ternal torque, b — 0 and choosing long enough ri so that in 
the long run Aflc relaxes back to zero while accommodating 
the AClc{t) > 0 values for the latest present observations. 
We have explored models with 0 ^ 6 1.1 x 10“^® rad s“^, 

corresponding to —7.2 x 10“^ ^ Next/N ^ 0. The values 
of b between 0.8 x 10“'^® rad s“^ and 1.1 x 10“^® rad s“^ 
yield reasonable fit results. Here we choose to explore the 
possibility of a permanent change in the external torque as 
reflected by 6 = 1 x 10“^® rad s“^. However, at present we 
cannot rule out fe = 0, a full decay. Parameters of the best 
fits with b — 0 and fe = 1 x 10“^® rad s“^ are shown in 
Table 1. The long term data display quasi-periodic resid uals 
with a period of ~ 400 days (lAntonopoulou et al.ll2014[) : we 
find a best fitting sinusoidal period P = 394 d by fitting the 
data from the last ~1500 days with a model involving only 
the terms that are dominant in the long term: the contribu¬ 
tions from inward moving vortices, the long-term offset with 
6 = 1 x lO”*^® rad s“^ and the sinusoidal term. In our further 
investigations comprising all the data we fixed this period 
for the sinusoid. The residuals of the 6 = 1 x 10“rad 
s“^ model also show initial fluctuations, which may be due 
to transient emission patterns in the magnetosphere. Future 
timing data will distinguish between these alternatives. 

To apply our extended creep model to the peculiar 
glitch of PSR J1119-6127, we use the spin-down rate data 
for the 2007 glitch, a total of 85 data points. The arrival 
time data from MJD 54268 indicate that a glitch has taken 
place since the previous data set on MJD 54220. The first 
post-glitch data fit to produce frequency derivative values is 
dated MJD 54300 dWeltevrede, Johnston fc Espinozall201ll l . 
The presently available spin-down rate and frequency data 
extending to MJD 56751 was kindly shared with us by 
Patrick Weltevrede ( P. We ltevrede private communication, 
lAntonopoulou et al.l (l2014fH . The time interval A between 
the actual glitch date and the first post-glitch frequency 
derivative values thus lies between A = 32 days and A = 80 
days. The coefficient I3/I ot the exponentially relaxing term 
is sensitive to the choice of A. We arbitrarily take A = 60 
days, which gives I3/I = 1.74 x lO”*^. 

We use the Levendberg-Marquardt method to And the 
best fit values of the paramete rs, starting from in itial guesses 
with MPFITFUN procedure (iMarkwardtl 1200911^ . The best 
fit is displayed in Figure 1 and its parameters are listed in 
Table 1. Inferred model parameter values corresponding to 
Eq. (15) are shown in Table 2. 

The long term remnant AHp of the glitch in frequency 
can be found by comparing the indefinite integral of the 
model for spin-down rate with the observed frequency resid¬ 
ual at the latest available data points. Using the frequency 
residual data on MJD 56688, AHp = 9.4 x 10“® rad s”*^ 
is obtained. We find AHp > 0, unlike the earlier result 
of Weltevrede et al. (2011), who found a negative long 
term frequency residual based on the latest post-glitch data 
then availabl e, but in agreement with th eir current estimate 
(Model A in lAntonopoulou et al.l (l2014fl '). 


^ http://purl.com/net/mpfit 
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Figure 1. Top panel: Fit to the post-glitch spin-down rate data 
with the model of Eq. (15), with A = 60 days and b = 1.0 X lO”*^® 
rad s“^. Middle panel: Zoomed version of top panel with model 
components representing contribution of exponential relaxation 
term (purple solid line), inward moving vortices (gray dashed line) 
and outward moving vortices (blue dash dotted line) are shown 
separately. Sinusoidal component and long-term offset b are not 
shown in the figure for clarity. Bottom panel: Difference between 
data and model. 


Table 1. Parameters of the best fits to the postglitch frequency 
derivative data following the 2007 glitch of PSR J1119-6127, with 
A = 60 days, and b = 1.0 X 10~^®raci s~^ (first column) and 
6 = 0 (second column) 


Parameter 

Value (Error) 

Value (Error) 

{ai)-i 3 {rad s~'^) 

1.68 (0.77) 

0.34 (0.29) 

{a2)-i3{rad s~'^) 

10.16 (1.43) 

5.93 (0.39) 

{a3)-\_2{rad s-^) 

9.92 (0.55) 

8.01 (0.38) 

Ctl 

-0.68 (0.15) 

-0.90 (0.07) 


6.34 (1.48) 

20.65 (3.82) 

Ti {days) 

1796 (211) 

20475 (14413) 

T 2 {days) 

159 (7) 

129 (4) 

T3{days) 

48 (2) 

58 (2) 

{b)-33{rad s~'^) 

1.0 

0.0 


6 DISCUSSION AND CONCLUSIONS 

We have examined the peculiar 2007 glitch of PSR J1119- 
6127 by extending the vortex creep model to take into ac¬ 
count (i) the possibility of a glitch associated change in the 
pulsar external torque, and (ii) inward motion of vortices. 
Both of these effects can be induced by a starquake that trig¬ 
gered the glitch. We model the peculiar glitch of PSR J1119- 
6127 as follows: a crustquake occurs, causing the crustal 
plates to move towards the rotation axis, together with some 


Table 2. Inferred Parameters with b = 1.0 X 10 ^^rad s ^ . 


Parameter 

Value 

(¥)_, 

1.11 

{¥)_. 

6.70 

{ 7 )_. 

1.74 

t'o{days) 

2046 

to {days) 

317 

(SQ ')_2 (rad s~^) 

2.69 

((5r2)_2 (rad s~^) 

4.16 

( AiVe„t \ 

V ^ 7-4 

-6.58 


pinned vortices. At the same time, some vortices affected by 
the crustquake are unpinned and move outward. The glitch 
is due to the angular momentum transfer associated with 
the sudden outward and inward vortex motions. Magnetic 
field lines, which move with the conducting crustal plate, 
change the external torque and give rise to the abnormal 
emission properties. 

In contrast to the changes in other pulsars’ glitches, 
the long-term change in spin-down rate, after transients 
are over, is (possibly) positive for PSR J1119-6127. In the 
creep process under the action of an external spin-down 
torque, the inward motion of vortices is thermodynamically 
unlikely, unless induced by a driving force such as arising 
from crustquake induced motion of crustal plates that carry¬ 
ing pinned vortices inwards. Inward vortex motion increases 
the lag between local superfluid and normal matter rota¬ 
tion rates from the steady state value, thereby accelerating 
rather than cutting off the creep process. This in turn in¬ 
creases the rate of angular momentum transfer to the crust 
and thereby decreases the spin-down rate of the crust, pro¬ 
ducing a positive change in the observed crust spin-down 
rate. By contrast, in standard glitches vortices move out¬ 
ward, decreasing the lag and turning off or suppressing the 
creep process which transfers angular momentum from su¬ 
perfluid to the crust; this leads to a negative step in the 
spin-down rate of the crust. Glitch induced steps of either 
sign arising from the offset in the vortex creep process al¬ 
ways relax back to the pre-glitch spin-down rate as the creep 
process heals back to the steady state. The model we have 
fitted to the spin-down rate data after the 2007 glitch of PSR 
J1119-6127 includes the creep response to outward and in¬ 
ward vortex motion as well as a glitch associated change in 
the external (pulsar) torque. 

lAntonopoulou et al.l ll2014h fit their data set with two 
models, model A with a long term exponential relaxation, 
and model B with a negative frequency second derivative 
(Ai>p < 0 in their notation). The positive Aflc decays to¬ 
wards zero in both models, asymptotically in the case of 
Model A. The two models leave comparable residuals. Our 
investigation of crust breaking, giving a permanent change 
in the external torque and spin-down rate and causing in¬ 
ward vortex motion, is complementary to their work. Future 
timing observations will decide if the offset in the spin-down 
rate is really permanent or relaxing; and searches for ther¬ 
mal signals accompanying future glitches of PSR J1119-6127 
will distinguish between the different models. 
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We must use the long term remnant of the frequency 
glitch in Eq. (13) to constrain Ib and Ib', the moments 
of inertia of superfluid regions which transfer angular mo¬ 
mentum only at glitches, due to the outward and inward 
motion of unpinned vortices, respectively. Using the val¬ 
ues of Sfls and SQ'g from Table 2 leads to the constraint 
9.4 X 10”®/ = 4.2 X 10"®7b - 2.7 x This gives 

Ib/I > 2.2 x 10“^. Then the superfiuid creep regions with 
a total moment of inertia Is > Ia + Ib + I a' + Ib' + h > 
Ia + Ib + I A' + ^3 > 20.4 x 10“^/ is effected by the glitch 
event. The region with moment of inertia I 3 = 1.74 x 10“^/ 
comprises most of the moment of inertia in pinned super- 
fluid. 

_ Recent calculations dChamell l2013l : lAndersson et al.] 

I2OI2II show that Bragg scattering of conduction neutrons 
from nuclei in the neutron star crust induces a neutron ef¬ 
fective mass that is larger than the bare mass. This “en¬ 
trainment” of superfiuid neutrons in the crust by the crys¬ 
tal lattice requires that the actual moment of inertia as¬ 
sociated with the superfiuid response is larger by a fac¬ 
tor m*/m„ > 1 where m* and m„ are effective and 
bare neutron masses in the lattice. The moment of inertia 
la ^ {rnnlm-n) 20.4 x 10“^/ associated with creep cannot 
be accommodated by the crust superfiuid alone for most 
neutron star models, even without the effective mass cor¬ 
rection. In addition to the crust superfiuid, other locations 
are required to sustain vortex creep. Contribution from vor¬ 
tex line-toroidal flux line pinning and creep a t the outer- 
core of the neutron star dSiderv fc Alpaj|2009lf . which has 
a comparable or larger moment of inertia than that of the 
crust su perfiuid, could provide the re quired extra moment of 
inertia dGiigercinoglu fc Alparll2014ll . The moment of iner¬ 
tia of the creep region where vortex motion is controlled 
by the toroidal arrangement of flux lines can amount to 
Itor/I ~ 2 X depending on the radial extension of 

the toroidal field in the outer core. Creep here is in the 
non-linear regime. As no glitch associated vortex motion is 
expected, the response to a glitch is exponential relaxation 
dGiigercinoglu fc AlDar|[2014l i . The relaxation time rtor — 50 
days for PSR J1119-6117 parameters is in line with our es¬ 
timate of T 3 = 48 days (and with r„i = 33 days for the ap¬ 
plication to the Vela pulsar , which has parameters simil ar 
to those of PSR J1119-6117: iGugercinoglu fc Alnarl d2014h l. 
So, we argue that I 3 /I = 1.74 x 10“^ reflects the moment 
of inertia associated with the toroidal flux line region of the 
outer core. The moment of inertia of the crustal superfiuid 
participating in the glitch, when the crustal entrainment cor¬ 
rection is included, is Is,crust > {m 7 lmn){lA + 7b + Ia') ^ 
(m^/mn) 2.98 X 10“^7. This total moment of inertia fraction 
can be accommodated in the crust in neutron star models 
with hard equations of state, if the mean value of m*/m„ 
to represent the crust s uperfiuid is not much larger than 
1. Th e recent work of IPiekarewicz. Fattovev fc Horowitj 
d2014ri shows that the neutron star crust may maintain 
larger moment of inertia so that the above constraint is eas¬ 
ier to be satisfied. 

The total numbers of vortices displaced in this glitch 
are determined by the superfiuid angular velocity changes 
SQs and SQ's using Eq. (9). The number of vortices that 
have moved outward is found to be is SNout ~ 1.3 x 10^®, 
while the corresponding number for inward moving vortices 
is SNin ~ 8.4 X These numbers are typical of all small 


or large glitches, from Grab, Vela and other pulsars analyzed 
so far in terms of vortex unpinning, indicating a particular 
scale of the glitch trigger. 

The glitch associated change in the external torque con¬ 
tributes a constant offset from the pre-glitch behaviour that 
remains in the spin-down rate after all post-glitch relaxation 
is over. This term, denoted b in Eq. (15), indicates a change 
in the external torque, which leads to a change in the spin- 
down rate through ANext/Next = ACi/Q+AI/I. The actual 
fractional change in the moment of inertia associated with 
a possible quake must be less than the observed glitch mag¬ 
nitude, so |A7/7| < Afl/fl ~ 10“® << |AS7/f7| ~ 10“'*. 
The measured permanent term b in Afl/fl therefore gives 
the fractional change in the external torque. Taking the ex¬ 
ternal torque to be essentially the dipole radiation torque, 
we have: 


ANsxt ^ A^ ^ AB± _ ^ AB± 
Next flc B± B± 


(16) 


as the term 3Aflc/llt: ~ 10“® is again negligible. We as¬ 
sume that the magnetic field change is associated with crust 
breaking, involving broken plates of size D distributed in a 
ring of the crust of width D and radius R cos a from the 
rotation axis, with each plate moving a distance D at the 
quake. The field moves with each broken piece of the con¬ 
ducting crust, without any change in the local field magni¬ 
tude, and orientation, which we take to be normal to the 
crust plate. The local field strength varies azimuthally in 
the broken ring. We further assume that the broken ring is 
in the polar regions of the magnetic field, so that the crust 
breaking has a strong effect on the external torque. This 
assumption is plausible if magnetic stresses play a role in 
crust breaking JFranco. Link fc EDsteirJl2nnnl : iLander et al.l 
l2014h . A schematic view of our model for external torque 
variation via change of magnetic field’s perpendicular com¬ 
ponent is depicted in Figure 2. The external torque variation 
is related to the change Aa in the angle between the rotation 
and magnetic axes: 


AB± Aa 
7?x tan a 


(17) 


From our estimate of the change in external torque, 
ANext/Nsxt — —6.58 X 10“* given in Table 2, we obtain: 


Aa = tana = (—3.3 x 10 *) tana (18) 

2 Next 

which is between (—1.0 x 10“*) and (—2.8 x 10“*). To 
obtain this, we have used the range o f a considered by 
IWeltevrede, Johnston fc Espino^ ll201ll l. a ~ 17° — 30° 
corresponding to an emission height 500 km and a ~ 
30 ° — 40 ° corresponding to an emission height 1800 km. 
The tiny change Aa in inclination angle cannot be resolved 
as an observable glitch associated pulse shape change in the 
present radio timing data. Motion of crustal plates towards 
the pole (Aa < 0) results in a reduction in the moment of 
inertia of the solid, and therefore an increase in the spin- 
down rate. This is the signature of a crust quake in a spin¬ 
ning down pulsar, tending to make the shape more spherical. 
The reason for the external torque change is likely to be a 
starquake, inducing motion of crustal plates, and reducing 
B± , as the surface magnetic field moves with the conducting 
plates towards the rotation axis. 

The fractional change in moment of inertia due to the 
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motion of the crustal plates is AI/I ^ {m/M)Aa Aa ~ 
10“'*, where m is the total mass of the moving plates and 
M is the mass of the entire star. Observing the direct effect 
of this actual change in the crnstal moment of inertia as a 
glitch is impossible. The glitch magnitude Afl is due to am¬ 
plification by the vortex motion triggered by crust-breaking 
and the resulting angular momentum transfer from super¬ 
fluid to normal matter. We have assumed that some of the 
vortices pinned to the moving plates are initially carried in¬ 
ward with the plates. This is possible when the increase Si}' 
in superfiuid rotation rate, due to the inward motion of the 
pinned vortices on the time scale of crust breaking is not 
sufficient for the local lag to increase from the steady state 
value to th e critical value, Sil' < ujcr — oJao , for typical val¬ 
ues of LOao (lAloar. Cheng fc Pineslfi989ll . These vortices will 
bend and are likely to be strongly perturbed by the sudden 
inward motion and become unpinned. The unpinned vortices 
will then move downstream azimuthally with the superfluid 
flow, causing more vortices to unpin and scatter outward 
until they reach a new radial position where they join the 
background vortex flow and creep. The avalanche of unpin¬ 
ning takes place rapidly on the glitch “rise” timescale. Some 
number SNin of vortices associated with the moving plates 
end up in radial positions inward of their original position 
while a number SNout of vortices end up in radial positions 
further out compared to their original position. The moment 
of inertia of the superfiuid creep regions affected by the in¬ 
ward motion of the plates and net inward vortex motion is 
of order 


I A’/I 


A-KpsR^D sin a cos^ a 
(2/5)Mi?2 

{2D/R), 


15/2 sin a cos^ a (D/R) 
(19) 


assuming a uniform density neutron star, and adopting 
sin a cos^ a = 0.3 for the range of a = 17 ° —40° indicated by 
IWeltevrede. Johnston fc Espino^ ll201lll . Using the value of 
I A' /I from fits, we obtain the ring width D ~ 6i?6 m, where 
is the neutron star radius in units of 10® cm. The number 
of vortices pinned to each plate is 

90 

SNpiate ~ -- 3 X lO^ilRl ~ 5.5Ri x 10®. (20) 

The total number of vortices associated with broken plates 
with a net inward motion during the glitch, is SNin ~ 8.4 x 
10*® vortices, as we estimated above from the results of our 
fits for to, a parameter independent from Ia'/^ which we 
used to estimate the plate size D. The number of plates 
involved should be ~ SNin/SNpiate ~ 10“*, in agreement with 
the number of plates in the broken ring, 2'kR/D 10*. 
A comparable number of vortices SNout ~ 1.3 x 10*® end up 
moving outward through a superfiuid region of comparable 
moment of inertia, I a- We find here an indication that the 
common scale, ~ 10*®, of the number of vortices unpinned 
in all pulsar glitches may be associated with the number of 
vortices in the typical plate size D involved in a triggering 
crust quake, multiplied by the number of plates involved, 
~ 2'kR/D 10*. These scales rest on the single parameter, 
the plate size D which must be related to the physics of 
the crustal solid. This plate size D is of the same order of 
magnitude as the ‘mountain’ height ~ 1 m estimated for 
the Crab pulsar dChamel fc Haensell 1200^ 1. Note that the 
critical strain angle 9cr at which the crust lattice breaks is 


6 or D/h where h is the radial thickness of the broken 

crustal plates. Thus, 

\1 m J \100 m J 

compatible with the results of lHorowitz fc Kadaul l|2009h for 
the critical strain angle. 

It is interesting to compare the moment of inertia frac¬ 
tions in crust superfiuid regions through which the un¬ 
pinned vortices moved during the peculiar glitch of PSR 
J1119-6127, given in Table 2, w ith the correspond ing Crab 
ll Alpar et al.l7199^ and Vela llAlpar et al.l Il993ll values, 
(0.01 — 1.87) X 10“® and (2.3 —3.4) x 10“® respectively, as this 
gives a lower limit on the moment of inertia fraction in the 
crust, lea ding to constraints on the neutron star equation 
of st ate dDatta &: Alpail 1 19931 : iLink. Epstein fc Lattim^ 
I1999II . With its characteristic age of ~ 1625 years, PSR 
J1119-6127 is between the Crab and Vela pulsars in age, 
but its implied crustal superfiuid moment of inertia frac¬ 
tion 2.98 X 10“® is comparable to the values inferred for the 
Vela pulsar. The qualitative evolution of glitching behaviour 
from Crab-like to Vela-like was proposed to be due to the 
development of connections in a network of vortex creep re- 
gion s, so that the momen t of inertia involved increases with 
age dPines fc Alpar! [198511 . PSR J1119-6127 should already 
have a sufficiently well connected vortex creep network. Pre¬ 
sumably the high magnetic field and associated stresses in 
the crust of this pulsar lead to high crust breaking activity. 
While similar moments of inertia in vortex creep regions I a 
are inferred, and the long term fractional offsets in the spin- 
down rate are similar in absolute value, for the Crab case 
no change in electromagnetic signature is observed. 

Magnetic stresses will play a role comparable to that of 
rotation induced stresses in conventional starquake models 
if, roughly, 

(^^0 - n") ~ \n\tp, 

where Bo and flo denote reference values of B and Q frozen 
into the crust. Using the glitch interval fg ~ 3 yrs, typical 
for young pulsars, we find that magnetic stresses can play 
a role where B > 10*® G in the crust. The magnetic field 
can have poloidal and toroidal components whose geometry 
will determine where in the crust the local stresses reac h 
the critical values for crust breaking dLander et al.l[2014l l. 
The higher multiple components of the magnetic field and 
the geometry of the stress tensor when both magnetic and 
rotational effects are included further complicate the situa¬ 
tion. In addition, the changes in electromagnetic signature, 
as seen only in PSR J1119-6127, are likely to occur if the bro¬ 
ken plate extends to the surface and leads to reconfiguration 
of the magnetosphere. This may explain why the behaviour 
exhibited by PSR J1119-6127 is rare. As a rough guideline, 
such behaviour may be exhibited by young pulsars with high 
magnetic field in young pulsars with high magnetic field. 

We should note the different responses of the crust and 
the superfiuid to a starquake. After the crust breaks and 
plates move towards the rotation axis in order to relieve 
their stresses, those broken pieces of the crust are stuck to 
new metastable positions, and do not come back to their 
pre-glitch sites. Thus, crust breaking and crustal motion 
are irreversible. Starquake induced inward motion of vor- 
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Figure 2. Starquake Model in cross section. The dotted area 
represents new position of crustal plates (broken ring) after the 
starquake. 


tices leads to a local excess of vortex lines and thereby to 
faster rotation of the superfluid. Creep becomes more effi¬ 
cient, evolving back towards steady state with an enhanced 
vortex current, as described by the vortex creep response 
employed in our fits. In the end, the superfluid relaxes back 
to the pre-glitch dynamical steady state, and all the parts of 
crustal superfluid and the rest of the star spin-down at the 
same rate. Indeed, without structural changes, the response 
of all internal torques will always relax back to Afic = 0. 

The main reason of the switch in emission patterns to 
intermittent and RRAT behaviour lasting for about a hun¬ 
dred days after the event in 2007 is likely to be the effect 
of the quake on magnetic field lines which are anchored to 
the crust. If the magnetospheric field pattern could move 
rigidly, without distortion, together with the motion of the 
crustal plate, as a consequence of the shift by Aq ~ 10“"^, 
there would be no significant change in the emission pat¬ 
tern. However, when a crustal plate moves in a quake, the 
elastic response of the field lines, which twist and reconnect, 
can amplify a small shift in the crustal position into a com¬ 
plex and drastic change in the emission pattern, helped by 
rotation which results in twisting of ma gnetic field fines an - 
chored to the highly conducting crust ( Beloborodov! 1200^ . 
The distorted magnetospheric geometry will subsequently 
relax towards a quasi-stable configuration for the new posi¬ 
tion of the plate. The changes in the emission pattern are 
observed for about a hundred days following the glitch, af¬ 
ter which the pulsar returns to its pre-glitch pulse shape and 
emission pattern. Twisting of field lines and their subsequent 
relaxation will also introduce temporary fluctuations in ar¬ 
rival times. Our timing model fits indeed leave relatively 
large residuals for about a hundred days in the post-glitch 
data given in the bottom panel of Figure 1. 

Although the re cently discovered mag netar anti-glich 
from IE 2259+586 bv lArchibald et ^ (l2013l ) has also shown 
a negative spin jump and changing emission features, the sit¬ 
uation is very different from that of PSR J1119-6127. In the 
case of IE 2259+586, the magnitude of the change in spin- 


down rate, |An| ~ 2.8|n| is too large to be associated with 
the superfluid regions in the star. It is likely that only a large 
change in the external torque is involved, as suggested by the 
violent change in emission, so that this magn etar ‘antiglitch ’ 
must be ex ternal/magnetospheric in origin (iLvutikov l2013l : 
lTondl2014h . 


In summary, the peculiar glitch of PSR J1119-6127 of¬ 
fers an invaluable opportunity for the reexamination and 
extension of glitch models to account for anomalous glitch 
signatures and transient emission phenomena initiated by a 
quake leading to a change in the external torque and trigger¬ 
ing the response of the superfluid regions of the neutron star. 
Our model predicts that the change in the external torque 
is permanent. The coincidence of the numbers of vortices in¬ 
volved in the glitch with the numbers inferred in Crab and 
Vela pulsar glitches is highly suggestive, supporting the ex¬ 
planation in terms of a crust breaking event with a typical 
plate size, which may be a common, even universal, trigger 
for glitches. If future timing observations rule out a perma¬ 
nent change in the external torque, this coincidence would 
turn out to be spurious. The explanation of the post-glitch 
evolution of Afic(t) in terms of internal torques responding 
to glitch associated inward and outward vortex motion, and 
relaxing to Aflc = 0 is viable when there are no observed 
changes in the external torque. 
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